APCTP-Pre209-003, KIAS-P09016, YITP-09-32 



Dynamical Compactification and Inflation 
in Einstein- Yang-Mills Theory with Higher Derivative Coupling 

Hironobu Kihara^, Muneto Nitta"'", Misao Sasaki^'*, Chul-Moon Yoo*, and Ignacio Zaballa^ 

Korea Institute for Advanced Study 
207-43 Cheongnyangni 2-dong, Dongdaemun-gu, Seoul 130-722, Republic of Korea 
"^Department of Physics, Keio University, 
O ' Hiyoshi, Yokohama, Kanagawa 223-8521, Japan 

o ; 

CN ■ # Yukawa Institute for Theoretical Physics 

^ ! Kitashirakawa Oiwake-Cho, 606-8502 Kyoto, Japan and 

I * Asia Pacific Center for Theoretical Physics 

r-| ' Hogil Kim Memorial Building # 501, 

p • POSTECH, San 31, Hyoja-dong, Namgu, 

^ , Pohang, Gyeongbuk 790-784, Republic of Korea 

^ ; (Dated: July 3, 2009) 
> ■ 

m ; Abstract 

4 °' Einstein-Yang-Mills theory in ten dimensions with a qnartic term 

' in the Yang-Mills field strength. We obtain analytically a class of cosmological solutions in which 

o : 

On . the extra dimensions are static and the scale factor of the four-dimensional Friedmann-Lemaitre- 

o . 

^ I Robertson- Walker metric is an exponential function of time. This means that the model can 

^ ' explain inflation. Then we look for solutions that describe dynamical compactification of the extra 

dimensions. The effective cosmological constant Ai in the four-dimensional universe is determined 
from the gravitational coupling, ten-dimensional cosmological constant, gauge coupling and higher 
derivative coupling. By numerical integration, the solution with Ai = is found to behave as 
a matter-dominated universe which asymptotically approaches flat space-time, while the solution 
with a non-vanishing Ai approaches de Sitter space-time in the asymptotic future. 
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I. INTRODUCTION 

There have been many attempts to consider extra dimensions in addition to our world 
of four-dimensional space-time, even though they have not been observed. The original 
idea dates back to Nordstrom [1] , Weyl [2] , Kaluza [3] and Klein [4] , who considered extra 
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dimensions in order to unify gravity and electromagnetic force in five space-time dimensions. 
Now the most promising unified theory, describing all fundamental forces including two 
types of nuclear forces, is considered to exist in ten, eleven or twelve dimensions after string 
theory, M-theory [5] and F-theory [6] have appeared. Superstring theory is consistent in 
ten-dimensional space-time. The extra six dimensions should be compactified. Some people 
require supersymmetry in four-dimensional space-time and the extra-dimensional space was 
assumed as a Calabi-Yau manifold. After the discovery of D-branes, D-branes or more 
generally "branes" offer the possibility of large extra dimensions or the brane- world scenario 
[7]. 

There were many efforts to describe cosmological solutions in the framework of higher- 
dimensional theories. Especially, the realization of de-Sitter like expansion of a 4-dimensional 
part has attracted much attention in connection with the inflationary scenario or the current 
accelerated expansion of our universe. One of those attempts is the flux compactification, 
which have received a lot of attention in recent years [8]. One of the most important and 
basic features of the flux compactification is to stabilize the size of a compactified space by 
certain configurations of high-rank differential form fields. 

Before string theory was discovered, Cremmcr and Schcrk studied an attractive possibility 
of compactification with the size of a compactified space being stabilized [9]. In order to 
achieve it they placed a non-trivial topological solution (soliton) of a gauge field on the 
compactified space, for instance a monopole on the sphere or a Yang-Mills instanton on 
the four-dimensional sphere S"^. In these cases, the compactified space is stabilized at a 
finite radius rather than decompactified to an infinite radius. So they called it "spontaneous 
compactification" . 

In this paper we would like to study if such a compactification can occur dynamically 
or not. In general, in order to stabilize a topological configuration of a Yang-Mills field in 
dimensions greater than four, we need higher order terms of the gauge field strength [10]. 
Some years ago Tchrakian introduced such a term, which we call the Tchrakian term, in order 
to generalize 't Hooft-Polyakov monopoles and Yang-Mills instantons to those analogues in 
dimensions greater than four [11]. The term is not renormalizable, but still quadratic in the 
time derivative. Recently some of the present authors have numerically studied a monopole- 
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like solution in six-dimensional Minkowski space by adding the Tchrakian term [12].^ One of 
the authors has further studied asymptotic solution of five-dimensional Tchrakian monopole, 
the generalization of Tchrakian monopole [14] . In the case of a six-dimensional sphere, an 
exact solution to a generalized self-duality relation has been constructed for SO (6) Yang- 
Mills fields with the Tchrakian term [15].^ Then this relation has been successfully embedded 
in the Einstein- Yang-Mills theory with the Tchrakian term in the geometry of the direct 
product of the four-dimensional Minkowski space (anti-de Sitter space AdS4) and of 
a constant radius, with (without) a ten-dimensional cosmological constant [17]. In this 
solution the gauge field distributes on homogeneously, so it is a natural generalization 
of Cremmer and Scherk [9]. At least for the Yang-Mills part, the configuration attains the 
minimum of the Bogomol'nyi bound when the radius of satisfies a certain relation with 
the gauge coupling constant and the coupling strength of the Tchrakian term. Therefore we 
expect that if we turn on the time variation of the space-time, we obtain a solution which 
describes the process of dynamical compactification. 

In this paper we consider cosmological solutions with a time-dependent scale factor of 
the three dimensions as well as with a time-dependent radius of S^, and study if there exist 
solutions with the radius of tending to a finite value, as a possible model of dynamical 
compactification . 

This paper is organized as follows. In Sec. II, we describe our theory, that is, the Einstein- 
Yang-Mills theory with the Tchrakian term in ten dimensions. We review the discussion of 
Bogomol'nyi completion. In Sec. Ill, we introduce an ansatz on the ten-dimensional metric, 
namely, the direct product of the four-dimensional Priedmann-Lemaitre-Robertson- Walker 
(FLRW) metric and with the radius as a function of time. Then we specify a gauge 
configuration which satisfies the self-duality relation and solves the Yang-Mills equation with 
the Tchrakian term. In Sec. IV, simple analytical solutions with a fixed radius of are 
given. The four-dimensional part of the solutions is either Minkowski or de Sitter, depending 
on the choice of the model parameters. In Sec. V, we consider solutions that describe the 
process of dynamical compactification. We investigate the behavior of the solutions both 
analytically and numerically. In general, the four- dimensional part behaves as de Sitter plus 

^ This was originally motivated by the computation of non-Abelian Berry's phases in T-dualized USp matrix 
model [13]. 

^ Generalization of instantons on the complex projective space CP^ has been also given [16]. 
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small oscillations, while the radius undergoes damped oscillations toward a finite value. 
For a particular choice of the model parameters that gives the product of a flat space-time 
times with a fixed radius, we find the four-dimensional part behaves as a dust-dominated 
universe, that is, with the scale factor proportional to t^/^. Sec. VI is devoted to conclusion 
and discussions. 



II. MODEL SETTING AND BOGOMOL'NYI EQUATION 



Let us start from the following action in ten-dimensional space-time: 



'S'tot ■— 'S'eh + ^YMT ' 
1 

IQttG 



^YMT--^ J {-FA *F + a\F A F) A *{F A F)-Vodv} . (1) 

Here dv is the invariant volume form, TZ is the scalar curvature with respect to the metric 
Qmn and F is the field strength two-form which takes values in the Lie algebra so(6). The 
star (*) denotes the Hodge dual operator acting on differential forms in ten dimensions. Our 
notation is summarized in Appendix A. For more details, see [18]. 

We consider the case where the space-time is locally a product space of A4 and jV. A4 
is a four-dimensional curved space-time and jV is a compact space. Let us denote the total 
space T. Metric on this space is 

ds'^ = g^i,{x)dx^dx'' + gij{x, y)dy^dy^ = ds%i + ds]^ , 

/.,i/ = 0,l,2,3, /,J = 4,5,---,9. (2) 

For the case where the field strength has only components along the compact directions, 
we can manipulate the Yang- Mills action as [17] 

— I Tr|-FA*F + Q;^(FAF) A*(FAF)} 
16 Jr 

= ^ y dv^^^ J Tr (^F^ ia-f7 *6 {F A F)j A *6 (^F ^ ia-fr {F A F)j 

± /" dv^^^ [ Tr 2ia-f7F A F A F , (3) 

where *6 represents the Hodge dual operator along the compact direction J\f. We call this 
procedure Bogomol'nyi completion. The term Q := fj^Tt^yF^ is a surface term and it gives 



the bound on the energy density. Then the Bogomol'nyi equation is 



F ^ ia-f7 *q{F AF)^0 . (4) 

If either of these equations is satisfied, the energy attains the minimum given by Q irrespec- 
tive of the sign ±. 

Suppose that ^4^°^ is a solution of equation of motion and F^^^ is the corresponding field 
strength. We denote the fluctuations around this solution SA, A — A^^^ +6A. Let us expand 
the left hand side of Eq. (4) in terms of these fluctuations: 



F - ia^7 hFAF^Bo + Bi{SA) + 62(6 A) , (5) 



where 



Bi{SA) -.^DoSA - ia-fY *q {DoSA A Fq + Fq A DoSA) , 

B2iSA) ■=qSA ASA- ia-f7 *6 {qSA A SA A Fq + Fq A qSA A SA) 

-ia-f7*6[{DoSA + qSAA6A) A{DoSA + q6AA6A)] . (6) 

Here Bq is the zero-th order term with respect to 6 A. The term Bi{6A) is linear in SA. The 
remaining B2{SA) includes higher order terms. By substituting this to Eq. (3), we obtain 

J Tr {-F A *F + a\F A F) A *{F A F)} 

~h> J ^^^^^Tr{i3i(M) A HBiidA)} + 0{dA^) + J (total derivative) . (7) 

Here the term Bq A *eBi{6A) is a total derivative term because A*-"-* is a solution of the 
equation of motion. The term 2Bo A *qB2{SA) includes indefinite quadratic form of SA, 
which might yield a tachyonic mass term. When A^^^ is a solution of ;Bo = which is one 
of Eq. (4), no tachyonic mass term appears in gauge sector. We mention that this does not 
necessarily mean the stability of the system under the presence of fluctuations of both the 
metric and the gauge field. This is an issue to be studied in the future. 
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III. ANSATZ FOR THE METRIC AND GAUGE FIELDS 



In this section we consider time-dependent solutions in the sense of Preund [19]. Namely, 
the metric is assumed to be in the form, 

ds^ = (is4 + dsl ; 

_ _ .^2 , r 2 2,^1 M^P . 2 _ r 2 2fe (a\ 

ds,- dt+L,e (i^^|^|2/4)2 ' ds,-L,e ^^^^^^y^y, (8j 

where the coordinates a* = {a^,a^,a^) span the three-dimensional space and — 
{y\ y^) span \a\^ {a^ + {a^f + {a^f and \y\^ := {y^f + {y'f + ■■■ + {y^f. 

The parameter k is ±1 or 0. 0i and 02 are functions of time t. Lq is a constant with 
dimension of length. The radius of 5"^ is given by it! = Lqc'^'^. This type of metrics was 
considered in various contexts, for instance in [20, 21]. 

The SO (6) gauge field configuration, represented in terms of differential forms, is assumed 
to be in the form, 

^ -7a5y"^V''+S (9) 



AqLoe'f'^ 

where a,b = 1, 2, ■ ■ ■ ,6 are the indices of the Lie algebra of S0(6), ^ab '■= (1/2) [Ta, 7b] are the 
infinitesimal generators represented by spinor, and is the vielbein of the six dimensional 
metric ds^, 

q is the gauge coupling constant. In the configuration of the gauge field A, the internal 
indices a,b,--- and the spatial indices I, J, - ■ ■ are identified by an embedding of the spin 
connection of the six-dimensional sphere into the gauge group. 
The exterior derivatives of the vielbeins are 

2(1 + MV4)^ (1 + MV4) 

^^^y'V^^VU^^dtAV^. (11) 



2Loe'^2 

Then the Ricci tensor components are given by 



= -3(01 + 0?) -6(02 + 0^) 

T^ij = 9ij (^<Pi + '^J^^ + 3</'i + 60102^ 

7e/J - (^02 + 5-^2^ + 60^ + 30102^ . (12) 
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The scalar curvature is 



n=64>i + 1202 + 120^ + 4202 ^ 3g^^^^ ^ 6 ^^^_2^^ ^ ^^_2^^^ ^ ^^^^ 

-^0 



Thus the Einstein tensor components are given by 

3 
LI 



Gu = Hi + 150^ + 180102 + ^ {ne-^^' + 5e-2^^) , 





= ^201 + 602 + 302 + 2102 ^ ^^20102 + (Ave-2^1 + Ibe-^^') 
Qij = -gij (^301 + 502 + 602 + 1502 + 150102 + (3«:e-2^i + lOe-^'f")^ . (14) 



As for the gauge field, its field strength is given by 

This satisfies the following duality relations, 

*F = -i/377 dv^'^'^ A F A F , *(F A F) = dv^^^ A 77 F , (16) 

where 

P .= Slie^<t>. . (17) 
3 

The self duality relation Eq. (16) becomes the Bogomol'nyi equation (4) ii P — a. In this 
case, there are no tachyonic modes at least in gauge sector. This determines a particular 
radius R — L^. oi the extra dimensions in terms of the gauge coupling constants q and a, 

Lc-.^J—. (18) 
V Q 

Because 02 depends only on the time coordinate, the exterior derivative of (3d^'^^v vanishes, 

d(/3(Z(^)v) = . (19) 
This means that the configuration satisfies the equation of motion, 

F)(*F) -a2£)|*(irAF) AF + F A*(FAF)} =0 . (20) 



The energy momentum tensor of the gauge field is given by 

a' 
3! 



'^MN — ^Tr ( —FmpFn^ + -^HmpqsHn^^^^ — ^SmnX > 



X := ^Tr (-^FM^^F^^ + ^HmnpqH'^''^'^ + ■ (21) 



Here Hukl are the components of F A F introduced in Appendix A. For our gauge config- 
uration we have 

_ 1 _ 1 



where 



In this gauge configuration, the Einstein field equations are 

= 30? + 150^ + 180102 + ^ («:e-^^^ + Se-^-^^) , (24) 

= (^201 + 602 + 30? + 2102 ^ -^20102 + -1 {Ke-^'t'^ + ISe"'^^)^ , (25) 

SttG \ —e-^'t'^ ( 1 - ^^e-^^^\ - -V^ 

[ 8q'Lf q^Lf ) 2 "^^J 

= - (^301 + 502 + 60? + 1502 + 150102 + (3fi:e-2<^i + lOe-^"^^)^ , (26) 

where the first equation is a constraint on the field and its derivatives, the Hamiltonian 
constraint equation, determining the three-dimensional hypersurface in the four-dimensional 
phase space. We note that the kinetic term in the Hamiltonian constraint is quadratic 
in the field velocities, and it has one positive and one negative eigenvalues. The above 
system of differential equations is invariant under the time translation and time reversal 
transformation. If «; = 0, there is in addition an invariance under the shift of 0i. The 
time evolution of the fields 0i and 02 is determined by Eqs. (25) and (26), describing the 
trajectory on the three-dimensional hypersurface defined by the constraint equation. 

It is convenient to express the field equations in terms of a rescaled time coordinate 
r = t/Lo, and introduce the following dimensionless parameters, 

SttG 

^2 7-2 



(27) 



q;2 



c = AtiGVoLI, (29) 

where c is related to the ten-dimensional cosmological constant A by c = ALq. In what 
follow, the r-derivative of a function /i(t) will be denoted by h' . For Lq = Lc the parameter 
h is fixed to the value h — 1/9, leaving only two free parameters in the field equations. 
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By manipulating the field equations we can reduce them to the following convenient set 
of two differential equations: 



= + 5(0^f + 6((/);)((/>^), (30) 

V2^cP'^ + 6{^',y + 3{^[m), (31) 
where Vi and V2 are defined by 

V^{<|>^, <t>2) X (^2) - Se"^^^ - ^e"^^^ , (32) 

= -ae-^<^^ (1 + 96e-^^^) + - - Se'^'^^ - ne-^^- , 
8 3 

T/2(02) := §e-^^^ (5 + 636e-^'^^) + | - Se"^^^ . (33) 

We can solve Eq. (30) for (f)[ to obtain 

<P[ = -30^ ± + 4(0^)2 . (34) 

By using this equation, we can eliminate (f)[ from Eq. (31). Then the Einstein field equations 
are reduced to a system of coupled differential equations given by 



cf>[ + 3cP',-^V,+A{<P',y = 0, (35) 
<P2 - 3(02)' + 3<PWVi+^{<f>'2y -V2 = 0, (36) 

where we have chosen the positive value of the square root in Eq. (35). In the next section, 
we look for a solution in which the extra-dimensional part of the metric is static, that is, a 
solution with ^2 — 0- In this case, 0i grows with time for the above choice of the square root 
sign, ensuring that the four-dimensional part of the metric describes an expanding universe. 



IV. SOLUTIONS WITH STATIC EXTRA DIMENSIONS 

In this section we consider sohitions in which the metric of the extra dimensional space, 
5"^, is static, that is when 02 =constant. In this case, Eq. (35) becomes integrable with 
respect to 0i, and Eq. (36) becomes an algebraic equation for e~'^'^^ . We note that we do 
not require our solution to satisfy the Bogomol'nyi equation (4). Hence for those solutions 
whose extra-dimensional radius is different from given by Eq. (18), the absence of tachyon 
modes is not guaranteed. Therefore we simply assume that there is a sufficiently wide range 
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of parameters in which there appears no harmful tachyons. This issue is left for a future 
study. 

Below we first consider general solutions. As we will see shortly, there is a particular 
solution given by e'^'^ = \J —kt + C for k = —1,0. Since this is somehow special, we treat it 
separately. 



A. The general case 

Static solutions of Eq. (36) are determined by the roots of 1^2(02) = 0. Let us set 
Z :— ae~^'^2 ^J\[Q note that Z oc (Lqc'^^^-s^ where Loe'^^ jg the linear scale of the extra 
dimensions. The equation V2 — becomes 

f{Z) = + 5iy,Z^ - 32iy,Z + = , (37) 

5 

where 

"^^=636=6W -2:=ac=^^. (38) 

Note that ui and 1^2 are independent of Lq and 02- As demonstrated in Appendix B, the 
equation f{Z) — has one or two real solutions Zi and Z2 (we assume Zi > Z2) provided 
that i/i and 1^2 satisfy a certain inequality. 

Let us first consider the solution Zi. The relation between the original variables and Zi 
can be written as 

o , , o SttG 1 
Llexp{2<P2) = Ll:=-^-, 39 

q Zi 

where Li represents the size of the extra dimensions. Thus the size of the extra dimensions 
is completely fixed by the coupling constants. 

As easily seen, Eqs. (35) and (36) are invariant under the rescaling, 

Lq^CLq , exp(02) ^ C^^ exp(02) , exp(0i) ^ C"^ exp(0i). (40) 



Using this degree of freedom, we fix the length scale Lq to be the size of the extra dimensions 
Li, or equivalently, we set (f)2 — 4 
a must be a solution of Eq. (37): 



Li, or equivalently, we set 02 = 02^^ = for this solution. Then we have Zi — a. Therefore 



5a 

/(a) =0<^c = 20 -— (5 + 636). (41) 
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With this normahzation, we find a positive real for c > 0. The condition Z2 < Z\ and 
c > give the folfowing inequahties: 



32 



> a > 



16 



5 + 636 - - 5 + 1266 
In terms of 02 these two solutions are given by 



"yf :=-;.log(Zi/a) = 



(42) 



(43) 



Those points are critical points or equilibrium solutions of the differential equation (36). 
The value of 02^^ is depicted as a function of a for each value of h in Fig. 1. The discussions 
in the rest of this subsection is valid for both solutions. 
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FIG. 1: Plot of <\>f^ as a function of o for 6 = 0.5/9, 1/9, 2/9 and 5/9. Because 0^^^ > 0^'^ = 0, 



(2) 



only the part of the curves above the line 
given by Eq. (42). 



(2) 



is meaningful. The value of a is bounded as 



Now we turn to Eq. (35). Setting 02 = 0, we have 



where Aj [i = 1, 2) is defined by 

A? :=^7rGLgx(0^^^)-5e-^^^" 



(44) 



(45) 
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We assume is positive. > gives an additional condition on the parameter a, 



4 16 

> a> 



(46) 



1 + 186 - " - 5 + 1266 
If this inequahty is satisfied, > 0, because ^nGLlx{(p2) — Se^^'^^ is concave downward as 
a function of e"^"^^ and its derivative at e~'^'^'^ = 1 is negative. The allowed region of a and 
b given by Eq. (46) is depicted in Fig. 2. 




0.0 0.5 



FIG. 2: Allowed region of a and b. The filled region is the allowed region which is bounded by the 
lines A? = and c = 0. 



The equation (44) can be integrated to give 

e'^' = 4- (e^'" + ne-^^^) , (47) 
2Aj 

where the origin of the time coordinate has been chosen to make the expression simple. 
Four- dimensional parts of these solutions are the same as those of Ishihara [20] . 

For large r, the term proportional to k, can be neglected and the scale factor of the 
four- dimensional space-time approaches R{t) = Loe^*"^, which describes a universe with 
accelerated expansion. Thus, although we do not claim that our model can give a realistic 
model of the universe, depending on the value of the constant Aj, it can reproduce a period 
of inflation in the very early universe or the present universe dominated by a very small 
cosmological constant. 
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B. The case Aj = 



When Xi — the solution (47) is no longer valid as it is, and we need a special treatment. 
In this case Eq. (44) imphes that k, must be either —1 or 0. In either case, (e"^!)' = yZ—n, 
and the solution is 



where C is an integration constant. 

The four-dimensional part of the solution for k = is flat. It was obtained in [17], 
which is almost the same as the one obtained by Cremmer-Scherk [9] , but with the radius of 

and the value of ten-dimensional cosmological constant modified by the presence of the 
Tchrakian term. 

The solution for k — —1 is also flat. The four dimensional line element is 



This metric covers the inside of either the future light cone or the past light cone of the flat 



V. DYNAMICAL COMPACTIFICATION 

In this section, we switch on the time dependence of (f>2 in order to see if our model has 
the possibility to describe the process of dynamical compactification. For this purpose, we 
analyze the stability of the solution 02 = 02^'* (= 0) and 02 = 02^^ in the second order 
nonlinear differential equation (36) in the case of k — both analytically and numerically. 

We first analyze the stabihty of the critical points analytically. For this purpose, we 
linearize the system of differential equations (see e.g. [22]). We find, however, that for 
Ai = this method is not sufficient to establish the stability of the critical point 02 = 02 = 0. 
Therefore we will try a different approach in this case. 

We first consider the critical point (02,02) = (0,0), which is a stationary or equilibrium 
solution of the differential equation (36). Denoting X :— 02 and Y := 02 and keeping only 
terms hnear in X and Y, Eq. (36) is written as the following system of first order differential 





(48) 




(49) 



space-time. 
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equations: 




where the matrix A is given by 





dV2 
d(j>2 



fc=0 




AX, 



(50) 



A 



dV2 
d4>2 



4>2=o 



1 

-3Ai 



1 

-3Ai 



(51) 



with 



A? = 3 - ^a{l + 186) , 
:=^(6-a). 



(52) 
(53) 



Note that Ai in the above is equal to the one defined by Eq. (45) with the normahzation 
condition (41). The solution {X,Y) — (0,0) is asymptotically stable if both of the two 
eigenvalues of the matrix A, 



Ai± Wa? + 



AdVo 



9( 



<t>2=0 



(54) 



have negative real part. 

To analyze the stabihty of the second critical point, (02, ^2) — (02^\ 0); we simply replace 
02 by 02 — 02^^ when linearizing Eq. (36). Then the eigenvalues are 



T2,± = 



J2) 



(55) 



For Ai = the real part of the two eigenvalues is zero, and the linear system corresponds 
to the harmonic oscillator for > 0. In this case we can not apply Poincare-Lyapunov's 
theorem above and additional information is required to establish the character of the critical 
point for the full nonlinear equation. We therefore treat this case separately. 



A. The case Ai > 



As in Sec. IV, we are interested in the solutions with Ai real and positive. Then the real 
part of the eigenvalues Ti -t is negative if 

dV2 



d(f)2 



> 0. 



(56) 



15 



This condition coincides with the condition Zi < which is satisfied when the parameters 
satisfy Eq. (46). Thus the critical point (02, ^2) — (0; 0) is stable. 

The system shows two different kinds of behavior in the neighborhood of the critical point 
(02,02) ~ (0,0)- When — | — y^i) < 0, the system undergoes damped oscillations 
with the amphtude decreasing as ^r'^l'^^^'^ , Otherwise the system is over-damped, showing 
simple exponential damping toward the critical point. 

For the second critical point (02, 02) = (02^\ 0), it can be shown that 

dV2 



(2) 16 o"^ aZ 

'1 



where /(2') is the function introduced in Eq. (37) and Z2 is the solution of /(^) = 
corresponding to the second critical point. This inequality follows from the fact that dj jdZ 
is a monotonically increasing function of Z with the unique df /dZ — at Z — Zq and 
Z2 Zq, which is proved in Appendix B. If df /dZ{Z'^ < 0, there are two real eigenvalues 
with opposite signs. Hence the critical point is an unstable saddle-point. In the special case 
when df/dZ{Z2) = 0, the first and second critical points become degenerate, and 02 = 
becomes the only equilibrium solution of the system. Note that we have = 21/2X1 
this case. 

To confirm the above stability analysis, we have performed numerical integration of 
Eq. (36). Our numerical results indicate that the linear analysis around the first critical 
point is accurate. In Fig. 3, we show the phase-space orbits of the solutions of Eq. (36) with 
Ai > 0. In this case, the first critical point is stable and the other critical point along the 
02 = axis is an unstable saddle point. The location of the saddle point depends on the 
values of the parameters a and b as well, and it roughly defines an effective stability radius 
for orbits near the solution (0,0). The time evolution of 02 for an asymptotically stable 
solution is shown in Fig. 7, where e'^'^ oscillates with a decreasing amplitude until 02 reaches 
zero. 

We have also integrated Eq. (35) for 0i. The time evolution of the three-dimensional 
cosmic scale factor R oc e'^^ is shown in Fig. 5. Initially when the oscillatory energy of 02 is 
non-negligible, the scale factor behaves as the one in a matter-dominated universe, R{t) oc 
T^/^. For sufficiently large r, after the amplitude of 02 has decayed exponentially, the universe 
eventually enters a stage of accelerated expansion, R{t) oc e^'^, with the (dimensionless) 
Hubble parameter H — Xi. 
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FIG. o: Phase space diagram 

{<l>2,<l>2) for Af = 5/12 and = 25/4. These are equivalent to take a = 1 and b = 1/9. The 
figure shows the two critical points of the system. The point (0,0)is stable while the second critical point is a saddle point with 
unstable orbits to its right. 

B. The case Ai = 

In the case Ai = the real part of Ti^± is zero, rendering the hnear analysis insufficient 
to determine the stability of the solution. Therefore we have to take into account the second 
order terms. 

To second order in X and Y, Eq. (36) gives the equations, 

X' = Y, (58) 
Y' = -Q^X + G (X, Y) . (59) 

where G{X,Y) is a quadratic function given by 

G (X, Y) = (^9^]2 ~ y ) + ~ 3rV4i]2x2 + 4^2 . (eo) 

Let us solve Eqs. (58) and (59) perturbatively. We assume > 0. To ffist order in X 
and Y, the system describes a harmonic oscillator. Namely we have 

X(r) = rcos(^]r + ^), (61) 
Y{t) = -nr sin(fir + ^/j) (62) 
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as a solution of the first order equations, where r and ip are arbitrary constants. Then the 
orbits in phase-space are elhpses about the critical point (0,0). 

Now we consider the effect of the second order terms. Here we just apply the so-called 
Krylov-Bogoliubov method of averaging [22] to study the behavior of the solutions.^ 

First, we introduce varying constants in the harmonic oscillator solution as 

X{t) = r(r) cos (fir + ^/'(r)) , (63) 
r(r) = -fir(r) sin (fir + ^(r)). (64) 

Then the system of differential equations may be expressed as 

r' = fr{T,r,ilj) , (65) 
^' = ^(r,r,^), (66) 



^ Detailed calculation is shown in Appendix C. 
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FIG. 5: This figure shows the time evolution of the scale factor R{t) for D.^ = 25/4 and Ai = 0.001. 



where 



fr{T, r,ilj) = - — sin {Qt + ip) G {r cos {Qt + ip) , -Qr sin {Qt + ip)) , (67) 
/^(r, r,ilj) =- cos {Qt + ip) G {r cos {Qt + ip) , -Qr sin (fir + ip)) . (68) 

Note that the right-hand sides of Eqs. (65) and (66) are periodic in r with the period 27iQ^^. 
Then instead of these equations, applying the Krylov-Bogohubov method of averaging we 
consider the time-averaged equations: 



f' 



27r 



fr{s,r,ip)ds = — 3fir 



^2 



iP' = — f^{s,f,^)ds = 



for r and ip. The solution is given by 



1 



3fir + const. 
Ip = const.. 



(69) 
(70) 

(71) 
(72) 
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These give approximate behavior of r and at sufficiently large r. Prom Eq. (63), approxi- 
mations to 02 and 02 for large r are given by 

02(t) ~ ^cosQt, (73) 
(l>2ir)---^smnT. (74) 

Then, for large r, Eq. (35) gives 

,, sin fir 2 

0i + TT- (^5) 

r or 

We can now read off an approximate solution for the field 0i, 

01 ~^logT-Si(fiT) . (76) 

Thus the scale factor behaves as 

R{r) = Loe^iM - LoT^/^e-'''^''^^ . (77) 

Apart from the small oscillations, this describes a matter-dominated universe. 

In Figs. 6-7, we show numerical solutions of the full non-linear system for Ai = 0. The 
numerical results are in good agreement with our analytical estimations. The time evolution 
of e'^i in Fig. 7 clearly exhibits oscillations around its central value r^/^ as we have shown 
analytically. 



VI. CONCLUSION AND DISCUSSION 



In this article, we studied time-dependent solutions of the ten-dimensional Einstein- Yang- 
Mills theory with the Tchrakian term. We obtained a class of simple analytic solutions in 
which the extra dimensions are static and the scale factor of the four-dimensional Friedmann- 
Lemaitre-Robertson- Walker metric behaves exponentially in time with the rate of expansion 
given by constants denoted by Aj (i = 1,2). Thus our model admits solutions describing 
inflation. 

We then considered a possible dynamical compactification of the extra dimensions by 
allowing them to be time-dependent. In the case Ai > 0, we found solutions in which the 
scale factor of the extra dimensions undergoes damped oscillations and approaches a constant 
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FIG. 6: The time evolution of _ l)/Lo with time r = t/Lg. We have chosen fl = 25/4, corresponding to the choice of 
a = 1. The oscillations in the proximity of the equilibrium solution <j)2 = are rapidly damped out as 1/t. 

value, while the four- dimensional scale factor approaches e'*'^'^. In the case of Ai = 0, we 
found numerically that the scale factor behaves as a matter-dominated universe R oc r^/^. 

Our model includes four dimensionful constants {G,Vo,q,a). They define four typical 
length scales in our model. Or if we fix the Planck scale or the gravitational constant, G, 
we are left with three dimensionless parameters. In addition, if we require the Bogomol'nyi 
equation to be satisfied, the linear size of the extra dimensions is fixed to be = ■^/Sa/q', 
and there remains only two dimensionless parameters. 

As is shown in Sec. II, when the radius of the compact direction is equal to L^., there are 
no tachyonic mode in the gauge sector. However, for a set of model parameters that gives a 
radius substantially different from Lc, a tachyonic mode may appear. To investigate when 
a tachyon appears and how it affects our model is certainly an important issue. Also for a 
complete analysis, in addition to fluctuations of the gauge field, it is necessary to include 
fluctuations of the metric and cross terms between them. These are left for future work. 

We also note that all the discussions given in this paper applies equally to the gauge 
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FIG. 7: The time evolution of e'^i/Lg for A = 0. In this plot Q = 25/4, which corresponds to a = 1. The time-averaged 
scale factor a{t) oc< e**! > describes a matter-dominated universe. 

group SU(4) in place of S0(6), because the matrices ■jab are block diagonalizable. Namely, if 
we project those matrices on the four-dimensional eigenspace with respect to the eigenvalue 
+ 1 of 77, we obtain self-duality relation of SU(4) without 77. Thus all cosmological solutions 
obtained in this paper are also valid for models with SU(4) gauge theory. Furthermore, since 
SU(4) is a subgroup of SU(A^), our cosmological solutions can be embedded into the Einstein- 
Yang-Mills theory with the Tchrakian term with SU(A^) gauge group. Generalization to 
other gauge groups like Eg or SO(A^) with N >8 remains as a future issue [23]. 

Recently some of us (HK and MN) considered the Bogomol'nyi equation on CP" [16]. 
By using the gauge configuration on CP^, we expect that we will be able to obtain similar 
cosmological solutions for CP^ compactification instead of studied in this paper. Also, 
it is interesting to see if similar cosmological solutions can be obtained for other types of 
compactification such as the compactification in terms of the Casimir energy [24]. These are 
also issues to be investigated in the future. 
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APPENDIX A: NOTATION 



1. Definitions and Properties of Tensors 

Here we explain our notation. The Einstein tensor and the energy momentum tensor are 
defined as 



1 ^2 33 

QmN '■= T^MN — -gUN^^ , TmN '■= — ^ — - ^JmN ■ i-^^) 



In terms of these tensors the Einstein equation is 

Qmn = ^t^GTmn ■ (A2) 

The Einstein tensor is obtained by the differentiation of the Einstein-Hilbert action S'^jj with 
respect to the metric g^^. The corresponding Levi-Civita connection Vj^p is defined as 

^NP-^lg^'^{dNgQP + dpgQN-dQgNp) ■ (A3) 
The Riemannian curvature 'R-fJpQ is defined as 

'^NPQ ■— <-'P^ NQ ~ <-'Q^ NP + ^ PA^ NQ ~ ^ QA^ NP ■ \-^^) 

The Ricci tensor TZmn and scalar curvature TZ are 

T^MN '■— T^'mqn ' ^ '■— g^^T^MN ■ (A5) 
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2. Differential Forms 



The tangent vector space of a point is spanned by Om- The basis dx^ of the cotangent 
space is the dual vector, dx^(9jv) = . For vector space V the Grassmann algebra A*{V) 
is defined as T{V)/I where T{V) is the tensor algebra T{V) := e~o^®*' and I is the two- 
sided ideal generated hy v ® v^v G V . We can define a linear operation which is called the 
Hodge dual. Let us fix p and q := D — p. The Hodge dual operator * is defined as 

By using the Hodge dual operation the metric on the differential Suppose that a; is a p-form, 

CO :^ ^UM,-M,dX^'--''' . (A7) 

The inner product is given by (00,00) :— co A *u!. Let us show the metric in terms of the 
component, 

— '^Mi-Mp<^Ki-Kp€ Ni-NqCi^ 



- ~ g(^p\yqi'^Mi-MpUJKi-Kpe "Ni-N.e "dv 

= -^a;M,..M,a;i.,..i.,A^--^-^--^- , (A8) 

where the metric ^^^-^pMi-i^p jg defined as follows: 

1 ^ 

^K,-Kp,Mv-Mp -.^-^Y^ sign(a) JJ . (A9) 

Here 6p is the p-th symmetric group consisting of all permutations of p characters. Finally 

we obtain 

uA*u^--UM,...MpUj'''-^- . (AlO) 
p\ 

The minus sign is from the fact that the signature of the metric qmn is Lorentzian. 



3. Clifford algebra 

We will use the Clifford algebra with respect to the six-dimensional Euclidean metric in 
order to represent the algebra so (6). Indices a,h — 1,2,-- - ,6 refer to the inner space. The 

24 



Clifford algebra is generated by 7a which satisfy 

{7a,76} = 25„6 , Tab ^[7a,7&] ■ (All) 

These generators are represented as 8 x 8 matrices. 7a6 satisfy the commutation relation of 
the Lie algebra so (6). Anticommutation relation of jab is 

{lal>,lcd} = 2jabcd-^S^d] ■ (A12) 

Here jabcd is an antisymmetric product of four generators defined as 

^' ae&p 

The chirality operator 77 is defined as 

77 = -1717273747576 , 7? = 1> 7? = 77 ■ (A14) 

By using this matrix, jabcd is written as a sum of products of 77 and jab, 

labcd = -^ea6cde/777e/ ■ (Al5) 

4. Notation for Gauge Fields 

The degree of freedom of a gauge boson is represented by the Lie algebra- valued one-form 

A, 

A:=^<7abci^^ F^dA + qAAA. (A16) 

where F is the corresponding gauge field strength two-form and q is the gauge coupling 
constant. Let us rewrite the action in terms of the components, 

^Tr i-F A *F) = i-Tr (FmatF^^) dv 

1 F^r,F"''^''^^jab7cddv 



4-32 



1 F^j,F'^''''''dv8{5aad-5M5ac) 



4-32 



4-2 



^ F^j^F'^'^'^'^dv . (A17) 
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For notational simplicity, we introduce the composite four form operator H, 

H:=FAF^ ^HMNPQdX^''''^ . (A18) 



The energy momentum tensor is 



1 -ab r^abP ^ 



2 



MAT — 7:-t^MP-^N 



T^HmpqsHn^'^^ — 7; 9 m NX 



where 



= ^IV (^-FmpFjv^ + ^HmpqsHn'''^'^ - \gMNX , (A19) 

X:=Tr(^-lFMivF^^ + ^i/MiVPgi^^^^^ + K)) • (A20) 
APPENDIX B: 02 :=CONSTANT SOLUTIONS 

In this Appendix, we derive an inequahty which gives the condition for Eq. (37) to have 
two real solutions. 

Because (ff{Z)/dZ'^ = 12Z^ + lOi/i > for arbitrary Z, the polynomial f{Z) has a 
unique minimum. This means that the number of real solutions of f{Z) = 0, Eq. (37), is at 
most 2. Let the value of Z at the minimum be Zq. Then f{Zo) must be non-positive for a 
real solution to exist: 

f{Zo) = Z^ + hv^Zl - 32z/iZo + \v^V2 < . (Bl) 

o 



Also Zq must be the unique real solution of the equation, 

= AZ^ + lOuiZ - 32i/i = . (B2) 



dZ 

Because — 32i/i is negative, Zq must be positive. In fact, by using Cardano's formula, we 
obtain 

/ ; \ 1/3 



Zo = I 4z/i + Av^\l 1 + ^^^ij - i^^^i + Aui^l + ^^Jyi I , (B3) 

which is manifestly positive definite. 

Now using df{Zo)/dZ — 0, the condition (Bl) reduces to 

-Z'q - 24Zo < . (B4) 
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Thus when the couphngs (G, Vq, q, a) satisfy the condition (B4), there are one or two real 
solutions Zi and Z2, {Zi > Z-i). Because Zq is positive, and we have the relation ^1 > ^0 > 
Z2, Zx is always positive if it exists. When the equality in Eq. (B4) is satisfied, we have 

Zi = Z2{= Zq). 

We assume that the parameters satisfy Eq. (B4). Then for z/2 > or equivalently c > 0, 
we have /(O) > 0, hence both solutions are non-negative: Zi > Z2 > 0. The solutions are 
given by the Ferrari's formula. 



Z^ei^ + eiVD; D := -7 (lOi/i + m) + 16ei^ . 
z 4 ^yu 



(B5) 



Here ei and 62 are ±1, and u is 
where 

^ ^ 25i^ 321/21^1 ^ 64^^2 _ 210^2 _ 250i^ , . 

3 5' 3 ^ 33 ■ ^ ' 

If i^i and U2 satisfy Eq. (B4), 7^/4 — 77^/27 > 0. This means that u is positive. The equation 
must have only one or two real solutions. This implies ei — 1 because D < if ei = —1. 
Thus the two real solutions are 



^1 = ^ + J-\ (lOi^i + u) + 16^ , ^2 = ^ - J— {lOu, + u) + 16^ , (B8) 



where Zi> Z2. 

The above expressions for the solutions Zi and Z2 are quite complicated as they are. 
However, using the scaling freedom of Lq, it is possible to simplify the expressions. For 
this purpose, let us first recapitulate Eq. (39) where the length Li representing the linear 
extension of the extra dimensions was introduced, 

.2 SttG 1 

= . (B9) 

Then we set Lq — Li, which implies Zi — a. 

Also for Z2, we may also simplify the expression in terms of a, b, c with the normalization 
Lq = Li. In this case, since Z — Zi — a is a, solution of f{Z) = 0, we have Eq. (41), 

5a 

/(a) = 0^c = 20-— (5 + 636), (BIO) 

o 
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and f{Z) can be divided by {Z — a). The quotient is 
63ab 



{{Z/af + {Z/af) + (4 - I) {Z/a) - | = . (Bll) 
In order to use the Cardano's formula, let us change the equation into the normal form, 

(Z/a + l/3)2 + A(Z/a + l/3) + S = , (B12) 

where 

636 ' 3 -6306 ^ ' 

This equation has only one real positive solution Z2. Therefore the solution is 

Finally let us derive the bounds on the parameters a and b. We assume c > 0. From 
Eq. (BIO), this gives a bound on a and 6, 

32-a(5 + 636) > 0. (B15) 

In addition, since Eq. (Bll) has only one real positive solution Z2 which is equal to or 
smaller than Zi, the left-hand side of it is non-negative at Z = a, 
63a6 /. c\ c 1 



4 



+ (4-0 = -(a(5 + 1266)-16) >0. (B16) 



Therefore the conditions that Z2 < Z\ and c > yield the bounds on the parameters a and 
h as 

>»>^^. (BIT) 



5 + 636 - - 5 + 1266 

APPENDIX C: ASYMPTOTIC BEHAVIOR IN THE CASE OF Ai = 

Here we derive the asymptotic behavior of the solution of the system given by Eqs. (58) 
and (59). Equations (65) and (66) can be written as 

^Qcos(J^r + V) + (^j + 6Q)Qsin(fJr + V) = -3^^ - f 6^^ - — ) cos2(Qr + V),(C1) 
r \ 2 / 

^ = ^tan-^(QT + V). (C2) 
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By eliminating the term ■^'/r from these equations, we obtain 

2 -?>n + T , (C3) 



r' 



- = 3Qt 

r 



r 

where 

-n - — ) sin(nT + V) + 30cos(20t + 2^^) - ( tt^^ - ^ ) cos(3nT + 3^'). (C4) 

2, oil J \z Oil J 

We can integrate this to obtain an expression for r, 

J drJ^ . (C5) 

As for the angle ip, from Eqs. (C2) and (C3), it satisfies 

ij;' = -rQ cos{nT + ip) (s + 6 sm{QT + z/;) + ^6 - cos=^ ((^r + . (C6) 

As clear from this equation, tends to a constant for r — > 0, Then T will be a function 
oscillating around zero. This imphes that the integral of in Eq. (C3) cannot be large. 
Thus in the region where r is large enough, r damps out in time as l/r, 

^ ^ 3Qr- JdrJ^ ^ 3Q7 ' ^^^^ 
This is consistent with our anticipation that ip tends to a constant. Therefore ignoring an 
irrelevant integration constant, the asymptotic behaviors of 4>i and 4>2 at large r are given 

by 

1 2 
02 ~ 777— sin l]r , 01 ~ -logr - Si(fir) . (C8) 
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